Pearling Instabilities of Membrane Tubes with Anchored Polymers 
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We have studied the pearling instabihty induced on hollow tubular lipid vesicles by hydrophilic 
polymers with hydrophobic side groups along the backbone. The results show that the polymer 
concentration is coupled to local membrane curvature. The relaxation of a pearled tube is charac- 
terized by two different well-separated time scales, indicating two physical mechanisms. We present 
a model, which explains the observed phenomena and predicts polymer segregation according to 
local membrane curvature at late stages. 



PACS numbers: 87.16.Dg, 68.10.-m 

Single-component phospholipid membranes have been 
the focus of intense interest in recent years, as the sim- 
plest model system of biological membranes [Q. The 
latter are highly complex systems, comprised of a bi- 
layer consisting of many types of lipids as well as a 
mesh of macromolecules such as proteins and polysac- 
charides, which participate in a large variety of cellu- 
lar processes. A natural step to mimic this complex- 
ity in a simple model system is to study the associa- 
tion of polymeric molecules with self-assembled single- 
component phospholipid membranes 

Experiments with polymers, which associate with 
membranes by anchoring^ have revealed changes in the 
bending moduli of bilayers of single-tailed surfactants 
1^ , and striking morphological changes in vesicles |§,§,D . 
Anchoring occurs by the penetration of a number of hy- 
drophobic side-groups grafted along a hydrophilic back- 
bone into a bilayer. Hollow tubular vesicles incubated 
in a solution of anchoring polymers having a polysac- 
charide backbone develop a pearling instability, above 
a threshold polymer concentration This instability 
was not observed with purely hydrophilic polymers, and 
was effected by hydrophobic groups alone (without the 
backbone) , but only at concentrations five orders of mag- 
nitude higher than with anchoring polymers. 

In this experimental and theoretical study we investi- 
gate the mechanisms responsible for pearling in our sys- 
tem. It has been suggested that the induction of curva- 
ture by the anchors, which sink into the membrane to 
a depth of half a bilayer, may constitute a mechanism 
which drives the pearling instability 1^. We present a 
novel nonequilibrium experiment, which shows that two 
independent mechanisms - spontaneous curvature and 
area difference - contribute to the pearling phenomenon. 
We also present new experimental and theoretical find- 
ings showing inhomogeneous shapes at late stages. 

In our experiments, vesicles were made of stearoyl- 
oleoyl-phosphatidyl-choline (SOPC) with Cis alkyl 
chains. The polymer used was hydrophilic dcxtran with 



a molecular weight of 162, QQQg/mol, functionalized both 
with palmitoyl alkyl chains Cie and dodecanoic NBD 
chains as fluorescent markers. The anchors are dis- 
tributed statistically along the backbone, spaced four 
persistence lengths apart on average (1 alkyl chain per 25 
glucose units). 1 jil droplets of SOPC in a 7:1 chloroform- 
methanol solution {10 mg /ml) were placed on a glass cov- 
erslip forming small lipid patches. The sample was pre- 
hydrated for 20 minutes under water-saturated nitrogen, 
then hydrated with 0.1 mM potassium buffer at pH 6.5. 
A large number of hollow tubes with one or more lamellae 
formed after hydration, connected to the lipid patch. The 
latter constitutes a reservoir with which the tubes can 
exchange lipid molecules. Experiments were conducted 
at room temperature, insuring that the membranes were 
in a fluid-like state. This allowed free diffusion of both 
lipids and anchored polymers along the bilayers. A drop 
of polymer solution of a given concentration was intro- 
duced through one of the cell sides. The polymer was 
added after the tubes were formed, and we therefore as- 
sume that it anchored mostly on the outermost leaflet 
of the membrane. Its concentration on the membrane 
grew from zero as more and more chains anchored from 
solution. Events were observed by phase contrast and flu- 
orescence microscopy and recorded on video. The NBD 
markers were excited with an argon laser (488 nm) , and 
observed with a CCD camera. 

Fig. 1 shows a tube undergoing pearling after addition 
of polymer. The instability typically starts near the end 
cap of the tube (Fig. la), and gradually propagates along 
the axis (Fig. lb). This is presumably since the polymer, 
diffusing from the side of the chamber, first reaches the 
tip of the tube. Existing pearls become gradually more 
spherical (Fig. Ic). The rate at which pearls form de- 
pends on the number of lamellae in the walls of the tube, 
as well as on the concentration of polymer in the wa- 
ter (2 — 100 mg/l), varying from several seconds to many 
minutes. A salient feature of Fig. Ic is the increase in 
pearl size towards the end of the tube. This size gradient 
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becomes much more pronounced at very long times, as 
Fig. Id illustrates for another vesicle. The string of pearls 
separates into a set of small, nearly-uniform spheres con- 
nected to a group of much larger spheres. 




FIG. 1. Snapshots of a multilamellar tubular vesicle un- 
dergoing a pearling instability: (a) 0, (b) 70 sec. and (c) 150 
sec. after onset of pearling. The concentration of polymer is 
below 2mg/l. (d) Image of an inhomogeneous pearled struc- 
ture at late stages of the instability, 900 sec. after onset of 
pearling. The polymer concentration is below IQQmg/l, but 
much larger than in (a-c). The scalebars represent 20 fim. 

We have studied the onset of pearling by measuring the 
wavelength of the fastest growing mode, P, just above 
threshold, as a function of the radius of the unpearled 
tube, i?o- Within experimental accuracy P is linear in 
i?o, with a slope k = 2ttRo/P = 1.07 ± 0.05. This re- 
sult agrees well both with the value k — \ predicted 
by theoretical models based on induced curvature as the 
driving mechanism , and with experimental results on 
pearling instabilities in single-tailed surfactant systems 
In tension-induced pearling k w 0.7 and thus 
our measurement excludes pearling due to polymer or 
flow-induced tension in our system. 

Polymers anchoring to one side of a membrane can 
induce curvature by two mechanisms. The first is an in- 
crease in the area of the outermost monolayer into which 
anchors sink . The second mechanism is a local defor- 
mation of the membrane which can be induced either by 
the anchors regarded as inclusions [jlTIl , or by an entropic 
pressure exerted by the polymer backbone [p^ . These 
mechanisms form the basis of two models that describe 
the tendency of a membrane to display curvature p3| : 
the area difference elasticity model (ADE) and the spon- 
taneous curvature model (SC), respectively. Calculations 
of equilibrium shapes of vesicles with cylindrical symme- 
try based on both the SC and ADE models indeed yield 
pearled shapes of constant mean curvature called Delau- 
nay surfaces [Q. The equilibrium shape is not sensitive 
to the pearling mechanism. 

Nevertheless, the dynamics of pearling may allow us to 
distinguish between the two mechanisms and reveal their 
presence, since they are characterized by well-separated 



relaxation time scales. Inhomogeneities in both area 
difference and spontaneous curvature decay diffusively 
jlSf . The relaxation of spontaneous curvature is asso- 
ciated with polymer diffusion in the membrane. The rel- 
evant diffusion constant was measured for various macro- 
molecules and falls in the range of 1 jirn^ / sec < Dgc < 
5 nw? /sec Area difference relaxes via the sliding of 
one monolayer with respect to the other, and does not 
involve diffusion of molecules over large distances. The 
diffusion constant associated with ADE can be estimated 
from dimensional analysis to be Dade ~ Ko/b, where Kq 
is the compression modulus of the membrane, and b is 
the friction coefficient between the two leaflets of a bi- 
layer. Experimental estimates of these parameters lead to 
50/iTO^/sec < Dade < 500 /im^/sec |jl^. Changes in the 
shape of a vesicle are coupled to motion of the surround- 
ing water, where the energy dissipated in the water is of 
the order of the curvature energy of the bilayer. From 
dimensional analysis we find that for both mechanisms, 
ADE and SC, the diffusion constant associated with the 
displacement of water is D^yd ^ k/tiRq ~ bOfim? / sec, 
where k is the bending modulus of the membrane and 
rj is the viscosity of water. Thus, it is difficult to dif- 
ferentiate between hydrodynamical and ADE effects by 
measuring the diffusion constant, but it should be easy 
to distinguish between the SC and ADE mechanisms. 

To study the relaxational dynamics, a micropipette 
was used to deliver locally a small volume (~ 10"'* /i^) 
of polymer solution of concentration 20 — lOOmg/Z close 
to a tube (Fig. 2). Fluorescence images show that the 
polymer was indeed concentrated in a local region shortly 
after the injection. Pearling occurred, after an induction 
time in which enhanced undulations were observed (Figs. 
2a-2d). The injection of the polymer was stopped after 
the pearled region increased to a length of order 100 /im. 
The pearled region then shrank gradually (Figs. 2e-2h) 
as pearls opened up one at a time starting with those far- 
thest from the region of polymer injection. This decay 
process took several minutes. Fluorescence images show 
that the amount of polymer in the surrounding water is 
negligible during the decay. 

We measured the length of the pearled region, L, as a 
function of time during the decay process. For diffusive 
decay, we expect k, 2D{t — t), where D is the rele- 
vant diffusion constant and i is a constant. Figure 3 is a 
typical example of the dependence of on time (mea- 
surements on several tubes yielded similar results). The 
system exhibits diffusive behavior at early times with a 
diffusion constant of Z) « 200 /im^/sec. We attribute 
this decay to a combination of ADE and hydrodynami- 
cal effects. At later times there is a sharp crossover to 
a much slower diffusive behavior with D ^ 5 ^m'^ /sec. 
According to our estimates, this corresponds to polymer 
diffusion, and is associated with the SC mechanism. Our 
results thus provide clear evidence that both ADE and 
SC mechanisms influence the pearling instability. 
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FIG. 2. Snapshots of a local pearling experiment, (a-d) 
show the formation of pearls as polymer is added from the 
micropipette at left, (e-h) show the subsequent opening of 
these pearls as polymer diffuses down the stem. Times are 2, 
5, 8, 11, 42, 46, 54, and 65 seconds after the pipette is brought 
to the vicinity of the vesicle. The scalebar represents 10 yLm. 

We now carry out a theoretical analysis of the pearling 
phenomena in the case of global application of the poly- 
mer (Fig. 1). We consider closed vesicles with poly- 
mer molecules only on the outer side of the vesicle, and 
assume that pearling is a result of the SC mechanism 
(All our predictions, except for the inhomogeneities in 
polymer concentration, apply equally well to the ADE 
model) . In contrast with standard curvature models [ p^ , 
the spontaneous curvature in our system is a local quan- 
tity. We assume that the spontaneous curvature is pro- 
portional to the polymer concentration on the membrane, 
pHq. p[r) is the fraction of the membrane area covered 
by polymer molecules at the position r, and takes values 
in the interval [0, 1]. Hq is the spontaneous curvature in- 
duced by full coverage of the polymer. We further assume 
that the observed vesicle shapes are close to the equi- 
librium shapes under the constraints of constant vesicle 
volume and membrane area. These shapes can be ob- 
tained by minimizing the free energy of the membrane 
in the presence of the polymer. The simplest free energy 
for our system is a sum of the curvature energy and the 
entropy of mixing of the polymer: 

F = j dA {2n{H ~ pHaf + 

^[plnp+(l-p)ln(l-p)]} , (1) 
where H is the local mean curvature, a is the character- 



istic linear size of an anchored polymer molecule and the 
integration is over the area of the membrane. In prin- 
ciple, the effects of gradients of polymer concentration 
should also be included in the free energy. However, the 
inhomogeneities at equilibrium do not induce an exten- 
sive free energy increase (see below), and we therefore 
ignore such terms. 
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FIG. 3. Squared length, , of the pearled portion of a 
tube as a function of time, t, after the pipette is brought 
to the vicinity of the vesicle (circles). The straight lines 
are fits of the early and late evolutions to the equation 
= 2D{t-t). At early times D « 200 ^m^/sec, while 
at late times D ~ 5 firn^ /sec. 

We have measured k w (20 ± 5) ksT using the pipette 
aspiration technique Structures with radii smaller 

than optical resolution ('^ 0.2 pm) were observed in our 
experiments, giving a lower bound of Hq > 10 pm^^ on 
the spontaneous curvature induced by the polymer. Fi- 
nally, a can be estimated as the radius of gyration of a 
polymer performing a two dimensional random walk on 
the membrane. Since the hydrophilic backbone is in wa- 
ter, we assume a good solvent in semi-dilute conditions, 
giving a radius between 40 and 80 nm. 

We consider very long, cylindrically symmetric vesi- 
cles and ignore the existence of the end caps, since the 
length of most of the experimental tubes is larger than 
their radii by two orders of magnitude. To find the equi- 
librium configuration of the system, the free energy (|l|) 
was minimized with respect to the vesicle shape as well 
as the local polymer density. This was done under the 
constraints of constant vesicle volume, membrane area 
and total number of polymer molecules. 

Ideally, the equilibrium configuration of the system 
would have a homogeneous polymer distribution and a 
curvature H = pHo everywhere. Such a configuration 
minimizes the energy and maximizes the entropy simul- 
taneously. This is indeed possible for a range of values 
of pHo. For example, a vesicle with volume-to-area ra- 
tio A can have the shape of a cylinder of radius 2A and 
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curvature H = 1/(4A). Another shape, having the same 
value of A, is a chain of identical spheres of radius 3 A 
connected by infinitesimally narrow necks. In this case 
H ~ 1/(3A) everywhere. In fact, the vesicle's curvature 
may take any intermediate value, 1/(4A) < H < 1/(3A), 
because for each of these curvatures there corresponds 
a Delaunay shape with the same value of A. On 
the other hand, it is not possible to construct a shape of 
constant curvature for H > 1/(3A) or for H < 1/(4A). 

In our experiment the polymer adsorbs onto the mem- 
brane gradually. Therefore, at the very early stages of 
the experiment p < l/(4A-ffo), and the shape of constant 
curvature which minimizes the energy is a cylinder of ra- 
dius 2 A and curvature H = 1/(4A). This is consistent 
with our experimental observations. In principle, we also 
have to consider possible inhomogeneities in membrane 
curvature. However, our calculations show that for small 
p such inhomogeneities only increase the free energy. 

At intermediate stages of the experiment the polymer 
concentration is in the range l/(4AiJo) < P < l/(3Ai?o)j 
and long vesicles are expected to have a Delaunay shape. 
Indeed, the shapes of vesicles we observe at intermediate 
stages of the experiment are similar to Delaunay shapes. 
Such pearled shapes have already been observed ||^ , and 
the importance of Delaunay shapes for this system has 
been discussed in the literature ]l^ , [T9t . 

The situation becomes more interesting at late stages 
of the experiment when p exceeds the value l/(3A_ffo)- In 
this case, the best Delaunay shape is a chain of identical 
spheres of radius 3A and curvature H = 1/(3A). How- 
ever, a detailed calculation (to be presented elsewhere) 
shows that a configuration consisting of a chain of small 
spheres connected by infinitesimal necks to each other 
and to one large sphere has the lowest free energy. 

The small spheres have polymer density pi and radius 
ri — \/{piHa). The energy of this subsystem vanishes, 
since its curvature is Hi = PiHq. The large sphere has 
radius r2 and polymer density p2, and plays the role of 
a reservoir for the excess volume and polymer molecules. 
The ratios r2 /ri and p2 / pi can be calculated numerically 
as functions of the average polymer density, the area of 
the system and its volume. For the values of the model 
parameters discussed above we obtained r2/ri > 10 and 
P2I Pi < 0.3. This ratio of the radii is consistent with 
the experimental configurations seen at long times, where 
chains of very small spheres coexist with few very large 
spheres, all connected by narrow necks (see Fig. Id). 

The strong inhomogeneity in polymer concentration 
predicted by the theory, is particularly interesting be- 
cause it may show a qualitative difference between pre- 
dictions of the local SC model and those of ADE mod- 
els. The latter does not differentiate between polymer 
molecules and lipids; i.e., exchanging polymers with a 
larger number of lipids of the same total area does not 
change the energy. One can always use such exchanges to 
turn an inhomogeneous polymer distribution into a ho- 



mogeneous one of higher entropy without changing the 
energy of the vesicle. Hence, ADE models predict a ho- 
mogeneous polymer distribution. In the local SC model, 
inhomogeneities in polymer distribution induce the same 
undesirable decrease in the entropy of mixing. However, 
inhomogeneities lower the energy, and our calculations 
show that for reasonable parameter values, the reduction 
in energy does lead to sizable inhomogeneities. We intend 
to measure the polymer concentration on the membrane 
to test this interesting theoretical prediction. 
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